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Abstract. Quantum spin-lattice systems in low dimensions exhibit a vari- 
ed ' ety of interesting zcro-tcmperature phases, some of which show non-classical 
ff . (i-G-, non-magnetic) long-range orders, such as dimer or trimer valence-bond 
order. These symmetry-breaking systems with localized valence bonds are 
referred to as valence-bond lattices (VBL) in this article. A review of our sys- 
Q ' tematic microscopic formalism based on a proper set of composite operators 
CJ . for the ground and excited states of the VBL systems is given. The one- 
dimensional (ID) spin-2 frustrated model is investigated in detail. Several 
possible VBL systems on the ID spin-1 chains, the 2D square and kagome 
lattices are also discussed. That our microscopic theory guarantees the ro- 
2 . tational symmetry of the VBL systems is emphasized. 



1. INTRODUCTION 

Spontaneous symmetry breaking (SSB) (the symmetries of a many-body Hamil- 
tonian not being preserved by its ground state) has always been a fascinating phe- 
nomenon in physics. Spin-lattice systems provide ample evidence for such SSB. Per- 
haps the most well-known example is the ferromagnetic (FM) Heisenberg models 
which have the classical ground state with all spins pointing in the same direction, 
say, along the z-axis, thereby breaking the rotational symmetry of its Hamiltonian. 

The antiferromagnetic (AFM) counterparts, however, prove much more compli- 
cated in quantum mechanics. This is well demonstrated by the fact that even on 
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a bipartite lattice the classical Neel state consisting of two alternating spin-up and 
spin-down sublattices is no longer the eigenstate of the Hamiltonian. Despite that, 
a number of AFM systems show a nonzero, albeit reduced, Nccl-like order. For ex- 
ample, the 2D spin-i AFM Heisenberg model on the square lattice, which has been 
under intensive study since the discovery of high-temperature superconductors, is now 
widely believed to possess in its ground state a Neel-like order which is reduced to 
about two-thirds of the classical value due to quantum correlations [1,2]. 

In addition to the FM phase and Neel-like AFM phase, quantum spin models 
also exhibit SSB of non-classical (i.e., non-magnetic) types in their ground states, 
such as the dimcrization in a spin chain with two adjacent atoms forming a spin- 
singlet valence-bond (VB), thereby breaking the chain lattice symmetry. Since the 
total spin vector of an isotropic spin system is also a good quantum number, those 
ground states with zero total spin vector certainly have no classical counterparts. Any 
long-range order in such a many-body ground state must be quantum mechanical in 
origin and the corresponding broken symmetry is likely to be the lattice symmetry 
because the rotational symmetry is preserved by such a ground state. 

Working in the sector of zero total spin vector has a long history. In fact, the 
seminal Bcthe ansatz, which provides exact solutions for the spin-i Heisenberg chain, 
was first proved by Hulthcn [3] through finite-size calculations within the frame- 
work of resonant valcncc-bonds (RVB), although the term "RVB" was not used then. 
Anderson extended this concept of RVB (originally due to Pauline) to frustrated 
spin-lattice systems [4], and later to the high-temperature superconductor materials 
[5]. More relevant to the present purposes, Majumdar and Ghosh [6] found that the 
perfect dimer VB configuration which breaks the lattice translational symmetry is 
the exact ground state of the important ID spin-^ Heisenberg chain at a particular 
ratio of nearest-neighbour and ncxt-nearest-neighbour coupling constants. In the last 
six years or so, the VB basis of low-dimensional quantum spin-lattice systems has 
attracted a lot of theoretical interest [7,8]. In particular, Affleck et al. [7] discov- 
ered that a homogeneous VB configuration is the exact ground state of a particular 
spin-1 Heisenberg-biquadratic chain. This finding sheds considerable light on the 
well-known Haldanc conjecture [9] on the nonzero gap in the spin-1 Heisenberg chain. 
(The term "valence-bond solid" was then used for such a VB state, although it has 
no conventional symmetry breaking [10].) 

An important rigorous result of quantum long-range orders is provided by a series 
of spin-s SU{n) {n — 2s + 1) chains (or the SU{2) spin-s chains with Hamiltonians 
which project out singlet states). This series of Hamiltonians has been solved by a 
mapping to the spin-2 XXZ chain which is integrable by Bethe ansatz [11]. By the 
same mapping, it has been shown that the ground state of the SU{n) model for any 
n > 2 breaks the lattice symmetry with a double degeneracy [12]. The exact values 
of the corresponding dimcrization order parameter have been obtained. In particular, 
the dimerization order parameters for the S'C/(3) and S'C/(4) models are reduced to 
about 42% and 68% respectively [12]. 

Although there are no exactly-known examples, it seems possible that sponta- 
neous trimerization, which is characterised by a sequence of spin-singlet states formed 
from three adjacent spins, may occur for some systems with integer spin quantum 
numbers since the trimer state is also a rather stable configuration. Recently this 
possibility has been discussed for the spin-1 Heisenberg-biquadratic chains over an 



extended region of the coupling constants [13]. Furthermore, the spontaneous dimcr- 
ization or trinierization may also occur in higher-order dimensionalities. For example, 
it has been proposed that the 2D spin-i frustrated Heisenberg model on the square 
lattice may show a column dimer VB order over a small but nonzero region of the 
coupling constants [14]. A more complicated dimerization picture was suggested for 
the spin-i Heisenberg model on the kagome lattice [15]. The trimerization of spin-1 
models on the kagome lattice is also a subject to be discussed in this article. 

For convenience, the term "valence-bond lattice" (VBL) is used in this article to 
represent collectively all those quantum spin-lattice systems in which the simple VB 
configurations (i.e., dimer or trimer, etc.) are localized with the broken lattice symme- 
try. One defines a perfect VBL as a regular array of isolated simple VB configurations 
on a lattice. In general, one expects that the perfect VBL is not the ground state of a 
given quantum spin-lattice Hamiltonian under consideration. But if the system pos- 
sesses a VBL long-range order in its ground state, the perfect VBL should be a good 
starting point. The quantum correlations can then be analysed on the basis of the 
perfect VBL. This same strategy was employed in 1952 by Anderson [1] in his AFM 
spin-wave theory, in which the quantum fiuctuations from the classical Neel state are 
described by collective motions of two sets of bosons. The VBL systems are not unlike 
the quantum Neel-likc systems, despite the complication of their ground states being 
in the sector of zero total spin vector rather than of zero total spin along the z-axis. 
Similar to the Neel-likc systems, where the quantum fiuctuations are described by 
the spin-fiip operators (i.e., spin raising and lowering operators) with respect to the 
Neel model state, a proper set of composite operators first developed by Parkinson 
in 1979 [16] is used to describe the quantum correlations with respect to the perfect 
VBL model state. A similar spin-wave theory can then be made by a bosonization 
scheme for those composite operators. But a more systematic approach is provided 
by a powerful microscopic quantum many-body theory, namely the coupled-cluster 
method [17], based on those composite operators themselves. The restriction in the 
sector of zero total spin vector is guaranteed by a very useful and important theorem 
of the CCM. 

Recently, Bishop, Parkinson and Xian [18] have successfully applied the CCM to 
a number of quantum spin systems, including the spin-i AFM Heisenberg model on 
the square lattice. In their analysis, the Neel state was taken as a model state for 
the anisotropic-Heisenberg AFM systems. Upon the Neel model state the many-spin 
correlations are incorporated via a so-called correlation operator consisting of the 
spin raising and lowering operators with respect to the model state. Within a well- 
defined systematic approximation scheme amenable to computer-algebraic techniques, 
they have obtained excellent results for the ground-state energy, excitation spectra, 
and staggered magnetization as functions of the anisotropy parameter. Their CCM 
analysis not only produces the numerical results which are among the best estimates 
available today, but also enables them to study the quantum phase transition of the 
anisotropic-Heisenberg systems in an extremely systematic fashion [19]. From these 
experiences, one expects that the CCM analysis should yield similar good quantitative 
results for the VBL systems. 

Because of its simplicity, our microscopic analysis for the spin-^ frustrated chains 
is first given. Then the same analysis is extended to other systems, including some spin 
models on the 2D square and kagome lattices. The outline of this article is as follows. 



Sec. 2 considers the few-body systems and introduces the corresponding composite 
operators and their boson transformations. Sec. 3 is devoted to the study of the ID 
spin-i frustrated model, firstly by the spin- wave approximation via a bosonization 
of those composite operators, secondly by the more systematic CCM analysis based 
on the composite operators themselves. Extensions of the same analysis to the other 
systems, including the spin-1 Heisenberg-biquadratic chains and the some 2D models 
on the square lattice and the kagome lattice, are discussed in Sec. 4. A general 
discussion is given in Sec. V to conclude this article. A brief proof of the important 
symmetry theorem of the CCM is given in the Appendix. 



2. FEW-ATOM SYSTEMS AND VALENCE-BONDS 

As outlined in Sec. 1, our microscopic theory for a VBL system is based on a 
set of composite operators which are defined according to the Hilbert space of the 
corresponding few-atom system. Our discussion here is restricted to the two-atom 
and three-atom systems. The boson transformation of those composite operators and 
the spin VB notations are also given. 

2.1 Two- Atom Systems 

A two-atom system, each with spin ^ , has four states, a singlet and triplet states. 
The singlet state can be written, in the obvious notation, as 

|0> = ^(IT,i>-U,T»; (2.1) 

and the triplet states are, respectively 

|1> = IT,T>, |2) = -^(|T,i) + U,T)), |3) = U,i). (2.2) 

In a matrix representation, one denotes each of these four states by a column matrix 
with a single nonzero element. Any operator in this Hilbert space can then be written 
as a 4 X 4 matrix. Following Parkinson [16], operator Amn is introduced as having 
only a single non-zero element in a (4 x 4) matrix, namely (?n'|Am„|n') = Smm'Snn'- 
All single spin operators can now be written in terms of these sixteen, namely 



-(^02+^20 + ^11-^33), 


4- 


= 2(^^02 - 


-A20 + A11-A33), 


(2.3a) 


-=(^30-^01+^21+^32), 


S2 


1 //. 


- A30+ -421+^32 ), 


(2.3&) 


-j={Ao3 - Aio + Ai2 + A23), 
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= V2^^^"- 


- A03+ A12+ A23). 


(2.3c) 



The inverse transformations are clearly nonlinear. Therefore Amn has been re- 
ferred to as a composite operator [13]. In particular, Aqq corresponds to the spin- 
singlet projection operator, namely Aqq = 1/4 — Si -82. 



We notice that Aiq (A30) is an operator which increases (decreases) s^^^^^ (= 
sf + S2) by one unit, while A20 leaves s^^^^^ unchanged; their Hermitian conjugates 
(i.e., transpose matrices) have the opposite effects. Since any of the triplet states can 
be generated by letting Ano {n — f , 2, 3) operate on the singlet state |0), Ano play the 
role of creation operators with respect to |0); their transpose matrices correspond to 
the destruction operators. Using the following algebra 

AmnAkl ^^ AmlOnkj (2-4) 

which follows by definition, it is easy to see that any Amn can be expressed by a 
product of Amo and Aon. 

For a two- atom system each with spin 1, the dimension of the Hilbert space is 
nine. The total number of composite operators Amn {in, n = 0, 1, 2, . . . , 8) is eighty- 
one. Among them are eight pairs of creation and destruction operators with respect 
to the singlet state. As in the spin-i case, if one chooses the the singlet state 

|0) = ^(|1,-1) + 1-1,1) -|0,0)), (2.5) 

the eight pairs of the creation and destruction operators with respect to |0) are then 
denoted as A„o and j4o„ {n = 1, 2, . . . , 8) respectively. 

2.2 A Three-Atom System 

Similar to the two-atom systems discussed above, in order to construct the com- 
posite operators of a three-atom system, one should list all of its states. Consider the 
spin-1 case. There are 3^ = 27 states. Ref. [13] lists all of them in detail. Again, the 
singlet state 

|0) = -^(|0,l,-l) + |l, -1,0) + 1-1,0, 1)-|0,-1,1)- 1-1, 1,0) -|1,0,-1)), (2.6) 

is the first state, and the nine states of stotai — 1 (stotai ^ Si + S2 + S3) follow, and so 
on, until the last state with all three down spins. 

As before, ^00 for the present case is also the spin singlet projection operator by 
definition, and can be written as 

^00 = -r5'l23(6 + 5123 - 'S'123), S'123 ^ Si • S2 + S2 • S3 + S3 • Si; (2.7) 

to 

the creation and destruction operators with respect to the singlet state |0) of Eq. (2.6) 
are similarly given by Ano and A^n {n = 1,2, . . . , 26). In Ref. 13, an approximation 
is made by truncating the Hilbert space from twenty-seven states to the first ten (i.e., 
restricting to the subspace of Stotai ~ and 1). In this subspace, there are only nine 
pairs of creation and destruction operators which can be easily managed. 

2.3 Bosonization of Composite Operators 

A bosonization scheme for a set of operators usually starts from a reference. The 
reference of the bosonization scheme (e.g., Holstcin-Primakoff transformation) in the 



conventional spin-wave theory [1] is either the spin-up state or spin-down state. The 
reference for our present purpose is clearly the singlet state |0) of the correspond- 
ing few-spin system. Therefore, the similarities between operators s^ and Aqo and 
between s+ (s^) and Ano (Aon) can be clearly seen. 

While s^- and s obey the usual SU{2) angular momentum algebras, from 
Eq. (2.4), it is easy to see that Amn obey the following pseudo-spin algebra, 

[Amn, Akl\ = A„iiOnk ^ AknOlva- (^-"j 

Therefore, Amn has also been referred to as a pseudo-spin operator. From Eq. (2.8) 
one can make the following Dyson-Maleev-like transformation, 

d-l 

^00 = 1 - 2_, '^n""! A„o = a^^OO, Aon = ««; Amn=am'^n, (2.9) 

n=l 

where m, n = 1, 2, . . . , d — 1 with d the Hilbert space dimensionality of the few-atom 
system, and where a„, a^ are [d— 1) sets of boson operators, obeying the usual boson 
commutation, [am, a^] ~ Smn- 

By definition, the singlet state |0) is the vacuum state of the bosons, namely, 
a„|0) = 0, n = 1, 2, . . . , d — 1. The physical states correspond to the vacuum state 
|0) and the {d — 1) states with only one boson excited. Furthermore, as the matrix 
elements between the physical and unphysical subspaces are equal to zero, the trans- 
formation given by Eqs. (2.9) is exact at zero temperature just as in the case of the 
conventional spin- wave theory [1]. 

2.4 Valence Bonds 

In the above analysis, the singlet state is always taken as the reference with 
respect to which creation and destruction operators are defined. This is the essence 
of our microscopic theory for the VBL systems. 

As is well known, spin singlet states can be conveniently represented in terms 
of VB which in turn can be expressed by Schwinger bosons [7,8]. Schwinger boson 
representation is given by the following transformation, 

s+:=^a+b, s-^ab+, s' = ^(a+a- b+b), (2.10) 

where a,a~^ and 6,6"'" obey the usual boson commutations. It should be emphasized 
that Schwinger bosons are used here purely for the notational purpose. They should 
not be confused with the bosonization scheme of Eq. (2.9). A spin-s state with 
s^ = m (— s < m < s) is written in the Schwinger representation as 

\s,m) = ^; ^ ^, ' > ^, \V), (2.11) 

where \V) is the vacuum state of the bosons. A spin VB between atoms i and j is 
defined by a number of the so-called VB operators [8] 

^J - ^tK 4^^ (2-12) 



Figure 1. Three VB configurations: (a) spin-i dimer, (b) spin-1 dinier and 
(c) spin-1 trimer. A single bond is defined by C^, of Eq. (2.12). 



acting on the vacuum state \V). Using Eq. (2.11), it is easy to see that the spin-singlet 
states of Eqs. (2.1), (2.5) and (2.6) can be conveniently written respectively, apart 
from the trivial normalization factors, as one-bond Cjt^jy), two-bond (C]t^)^|y), and 
three-bond Ci2C2zCti\^) configurations. 

A general VB configuration can be easily drawn for a spin-s many-spin system. 
In Fig. 1, two dimer and one trimer configurations are shown. A many-spin ground 
state with zero total spin vector [stotai (= X^i ^i) — 0] is in general given by a linear 
summation of all independent VB configurations in which each atom is linked by 
2s VBs [7,8]. In general, different many-spin VB states are not orthogonal to one 
another. This makes working in the VB basis very difficult. But the ground state of 
some interesting quantum systems is dominated by a particular VB state consisting of 
an array of independent simple VBs such as those shown in Fig. 1 . These arc the VBL 
systems defined in Sec. 1. In the following sections, a systematic microscopic theory 
is developed by taking these perfect VBLs as the reference state and by employing 
the creation and destruction operators Ano and Aqk with respect to this reference. 



3. THE SPIN-i FRUSTRATED CHAINS 
3.1 Spin- Wave Theory 

The ID spin-i frustrated model is perhaps the simplest model with spontaneous 
dimerization. The model consists of A^ atoms each with spin ^ on a chain with nearest- 
neighbour and ncxt-nearest-neighbour interactions. The Hamiltonian is simply 

N 

H = ^{si ■ Sj+i + Jsi ■ Sj+2), (3.1) 

i=l 

where J is the coupling constant, the usual periodic boundary condition is assumed, 
and even A^ and a unit lattice spacing are also chosen for convenience. At J = 0, H 
is the well-known Heisenberg model which was solved exactly by Bethe ansatz [3] ; its 
ground state is gapless and has no long-range order. At J = 1/2, the ground state is 
given by the dimer VB configuration as shown in Fig. 1 (a) with a double degeneracy 



[6], 

N/2 

\D) = n |0>2.-i,2., (3-2) 

where the notation |0)j represents the singlet state of the pair given by Eq. (2.1). 
Let r denote each dinier in Fig. 1 (a), and Si(r) and S2(r) the two spins of the dimcr, 
Eq. (3.1) becomes 

JV/2 

ff == ^ [si (r) • S2(r) + S2(r) • Si (r + 1) + Jsi(r) • Si (r + 1) + Js2(r) ■ s^ir + 1)] . (3.3) 

r=l 

One can then express H in terms of the composite operators Amn by Eqs. (2.3). 

As discussed in Sec. 2, since the fluctuations with respect to \D) can be described 
by operators AJ^q and A^^ (n — 1,2,3), one can derive the equations of motion for 
all of these three sets of pairs. By employing the usual decoupling approximations 
and taking Ago ~ 1, it is easy to derive the spin- wave spectra (i.e., eigen modes). 
Parkinson [16] employed this method to obtain the triplet spectrum for the Heisenberg 
model (J = 0). 

Application of bosonization scheme not only provides a more systematic means to 
obtain the excitation spectra, but also allows one to study the ground-state properties 
as well. By Eqs. (2.9), one can further express H in terms of the three sets of bosons 
(a polynomial up to sixth order) . Diagonalization of the quadratic parts of H by the 
usual Bogoliubov transformations, one can easily obtain the ground-state energy Eq 
and excitation spectra ujq within the spin-wave approximation. They are given by 
respectively 

(3.4) 



Eq 3 r dor I ^ T 3 

^-ii ^[v/l-(l-2J)cos2,-l]--, 



and 



= v/l-(l-2J)cos2g. (3.5) 



Eq. (3.5) agrees with that of Parkinson at J = [16]. A discussion of these results is 
left to the end of this section. 

3.2 The Coupled-Cluster Method 

The COM has been successfully applied to a wide range of quantum many-body 
problems in both physics and quantum chemistry [17]. The interested reader is re- 
ferred to Ref. [17] for the general formalism of the CCM and to Ref. [18] for its 
particular application to the spin systems with the Nccl-likc order. Here its extension 
to the VBL systems is considered. 

(a). The Ground State 

The CCM ansatz for the ground ket state is \9g) = t:'^]'!'), where |$) is the 
so-called model state which is usually chosen as an uncorrelated many-body wave- 
function, and where S is the many-body correlation operator consisting purely of the 
creation operators with respect to |$). For the VBL problem under consideration. 



it is quite natural to choose the perfect VBL state as the model state. The creation 
operators with respect to this model state are clearly given by any combination of 
those operators A^q with n = 1,2, . . . ,d — 1, where d is the dimensionality of the 
corresponding few-spin system and r denotes its vector position in the VBL. Their 
Hermitian conjugates Ag„ are the corresponding destruction operators. 

The Schrodinger equation of the ground state, after a simple manipulation, can 
then be written as 

e-^iJe^|$) =^g|$), (3.6) 

where Eg is the ground-state energy, and where the similarity-transformed Hamilto- 
nian can be expressed as a series of nested commutators, namely 

e'^Hc^ = H+[H,S] + ^^ [[H, 5] , 5] + • • • , (3.7) 

which usually terminates at the fourth-order for most Hamiltonians with pair- 
interaction potentials [17,18]. For the present case, the series of Eq. (3.7) terminates 
because Hamiltonians always contain a finite-order polynomial of the destruction op- 
erators. 

Now let us focus on the spin-2 dimerization. The model state is the given by the 
dimer state, |$) = \D). There are three sets of creation operators, namely AIq,A2q 
and ^30. If one restricts to the sector of zero s^^^^j (= ^^ sf), the correlation operator 

S = J2n=l^n, with 

N/2 A'/2 |. 

5i = 2_^^rA2o, S2 = 



r-r' 



c(l) A-r Ar' c(2) Ar Ar' 

'-'r,r'^10^30 9! »','• 20^20 



N/2 

r r' .r" 



(3.8) 



1 



c(l) Ar Ar' Ar" ^ c(2) Ar Ar' Ar" 



etc. In Eq. (3.8) the primes on the summations imply exclusion of the terms with any 
pair of indices being equal. 

The ground-state energy is obtained by taking the inner product of the 
Schrodinger equation (3.6) with the model state \D) itself, namely 

Eg = {D\c-^Hc^\D); (3.9) 

and the correlation coefficients {Sr.r'....} in Eqs. (3.8) are determined by the coupled 
set of equations obtained by taking inner products of Eq. (3.6) with states constructed 
from the corresponding destruction operators, namely 

{D\Al^2C^'^He'^\D) = 0, Vr, (3.10) 

for the one-body equation; and 

{D\Al,X3^-^Hc^\D)=0, (i?|ylS2<c-^i/e^|i?) = 0, W,r'{^r) (3.11) 

for the two-body equations. The three-body equations and higher-order many-body 
equations are obtained in a similar fashion. 



The exact energy equation (3.9) can be straightforwardly derived as 

^ = i[(l-2J)(25(^)+6f)-a)-3], (3.12) 

where the lattice symmetry is used to set accordingly Sr — a, Sri\r2 — SrJ^n = 
br , with 1 — 1,2 and r = r2 — ri . 

The exact one-body equation (3.10) can also be easily derived. It couples only 
to the two-body and three-body cocfhcicnts. Similarly, the two-body equations (3.11) 
couple only to the one-body, three-body and four-body coefficients, and so on. One 
clearly needs to employ an approximation scheme for a practical calculation. The 
most common approximation of the CCM is the SUBn scheme which retains up to 
n-body correlation operators. Here the SUB2 scheme is considered, namely S —^ 
S'suB2 = Si + S2 and Sn — for n > 3. The one-body equation (3.10) yields an 
interesting solution, a — 0, implying no one-body correlations for the dimerization 
problem. Furthermore, the two-body equations (3.11) provide a solution in which the 
two sets of two-body coefficients are identical, namely 

bi^^ ^bi^^=br. (3.13) 

We notice that the model state \D) is in the sector of Stotai = 0, and the one-body 
correlation operator 5*1 will take the state out of this sector. We also notice that the 
two-body correlation operator S2 commutes with Stotai if and only if Eq. (3.13) is 
satisfied. All these imply that the ground state in our SUB2 approximation remains 
in the sector of Stotai — despite the fact that we started with operators in the sector 
of Stotai ~ '-'• ^^ ^^*^^' ^^^^ ^^'^^ property also holds at higher-order approximations in 
the above CCM analysis. Appendix provides a brief proof for a general theorem which 
states that the CCM coupled equations [e.g., Eqs. (3.10) and (3.11)] at any level of 
approximations always provide at least a solution which guarantees the symmetry of 
the model state if this symmetry is one of those belonging to the model Hamiltonian. 
This is certainly a big advantage because it is much more difficult to work in the 
sector of Stotai = than of s^^^^^ = 0. 

The energy equation is now reduced to 

^ = ^[(l-2J)&i-l], (3.14) 

and, after simplification, the two equivalent two-body equations in the SUB2 scheme 
are given by 

i J2 (KsSrp + K^br - 2Kibr+p + i^l ^ K'br+p-r') =0, r ^ (3.15) 

p=±l r'=/^0 

with Xi = 1 - 2 J, K2 = A{l-2Kibi), &uA K3, = Ki{l + Abl)-2{l + 2J)bi. A simpler 
approximation can be made from Eq. (3.15), namely the SUB2-2 scheme which retains 
only the single coefficient, 5i. Eq. (3.15) then reduces to 

1 - 2J + 2(3 - 2J)bi - 9(1 - 2J)bl = 0, (3.16) 
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which is easily solved. The full SUB2 equation (3.15) can also be solved analytically 
by a Fourier transformation method in a similar fashion as described in Ref. 18. Here 
only the final result is given by the following self-consistency equation for 5i, 

^'" SiT (^"^2^ / dq^/l-kiCos2q + k2Cos-'2q\ , (3.17) 

where the constants fci and k2 are defined by 



^,{AK,K2 + %Klh,-AKlXl fc.^l^ll^^^) 



k^ = —^{AKiK2 + mih^-AKiX), k2 = i^-^2 i^, (3.18) 



and where X = X]r=i ^r&r+ij which can be calculated self-consistently as 6i of 
Eq. (3.17). After &i is determined as a function of J, the ground-state energy is 
obtained by Eq. (3.14). Again the discussion of these results is left to the end. 

(h). The Excited States 

The CCM ansatz for the excited state is |*e) = X\^ g) = Xe^|$), where |*g) 
is the ground state as determined above and X is the excitation correlation operator 
consisting only of the creation operators as S does. Using the Schrodinger equations 
for the ground and excited states, one obtains 

e-^[iJ,X]e^|$) =eA:|$), e = E^-Eg. (3.19) 

In the so-called SUB(2,1) scheme, one retains up to two-body correlations in S 
and one-body correlations in X, namely, 5* -^ 5*31162 = 5*1 + 5*2 and X ^ Xi. There- 
fore one writes, Xi — J^r'^^^in- '^^•^ other two one-body excitation operators are 
given by replacing A^q by ^20 ^^'^ ^30 respectively. The coefficient AV is determined 
by the inner product of Eq. (3.19) with the state ^5^q|$). A Fourier transformation 
readily yields the following excitation spectrum with a lattice momentum q, 

K2 I 

e^ = — — yl — fc2Cos2g + fc2Cos2 2(7, (3.20) 

where the constants K21 k\ and fc2 are as defined before. The other two excited states 
with operator Ajq and AgQ produce the same spectrum. Therefore Cg is a triplet 
spectrum as expected. 

3.3 Discussion 

Fig. 2 shows the results of the ground-state energy per atom as a function of J 
from the spin-wave approximation, and from the SUB2-2 and full SUB2 schemes of the 
CCM. The numerical results [20] of Tonegawa and Harada, obtained by extrapolating 
the finite-size calculations for J < 1/2, and the exact results by Parkinson of the 
iV = 20 system for J > 1/2, are also included for comparison. At J = 1/2 (i.e., the 
Majumdar-Ghosh point) both spin-wave theory and the CCM approximations give the 
exact result of —3/8. This is not surprising because both take the dimer state \D) as 
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Figure 2. Ground-state energy per atom as a function of J. Shown are 
ttie results from tlie spin- wave tlieory (dotted), the SUB2-2 scheme (short- 
dashed), and the full SUB2 scheme (long-dashed). The terminating points 
are indicated. The numerical results from Ref. [20] are also included (solid). 



their model state. At J = (the Heisenberg point), spin- wave theory yields —0.4498, 
while the SUB2-2 and full SUB2 schemes yield -0.4268, -0.4298 respectively. They 
all agree with the exact result of —0.4432 by the Bethe ansatz [3]. But we notice that 
at J = spin-wave theory produces divergent results for other physical quantities 
such as the dimerization parameter discussed later. Furthermore, as can be seen 
from Fig. 2, the energy curve of spin- wave theory is symmetric about J — 1/2, while 
the extremely simple SUB2 scheme gives much better results for a wide range of the 
coupling constant J when compared with the numerical results of Ref. [20] . 

We notice that both spin- wave theory and the SUB2 scheme have two terminating 
points Jc , beyond which, namely for J < Jc and J > Jc , there is not real 

solution. For spin- wave theory, the two points are given by Jc = and Jc — I, 
while Jc = —0.4443 and Jc = 1.591 from the SUB2 scheme. The corresponding 
energy values of the SUB2 scheme are —0.5172 and —0.6977 respectively. In the 
past we had identified the SUB2 terminating points as the phase transition critical 
points for the anisotropic-Heisenberg models [18]. This is strongly supported by 
the calculations of the spin correlation functions and order parameters within the 
same CCM analysis. The following discussion of the triplet spectra of the spin-wave 
excitations also supports that the two terminating points J = Jc of the present 

dimerization case may again correspond to the quantum phase transitions of the 
frustrated systems. 
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The triplet spectra of Eq. (3.5) and Eq. (3.20) have the same quahtative be- 
haviour. Fig. 3 shows the schematic plot of the spectrum from the SUB(2,1) scheme 
of the CCM at several values of J . The spectrum clearly shows a nonzero gap between 
the two terminating points and the gap collapses at both the terminating points. In 
particular, the triplet spectrum is flat with a gap value of 1 at J = 1/2. This flatness 
implies no coupling between pairs of spins (dimers) at J = 1/2 within the approx- 
imations. Since the simple dimer state \D) is the exact ground state at this point, 
the two-body correlation can be easily included in the excitation operator X and 
the corresponding two-body coefficients can be determined by a simple variational 
procedure similar to the well-known Feynman theory for the excitation spectrum of 
the *He superfluid [21]. Hence the excitation operator with a lattice momentum q is 
written as, 

X, = Y, c^^^^Al, + ^ /,(r, r')A^o<- (3.21) 

r r.r' 

Taking ,fq{r, r') as the variational parameter to optimise the expectation value of the 
Hamiltonian, it is found that the gap in the spectrum is reduced by half at g = 
and TT, but remains at 1 at g = 7r/2 [22]. (The preliminary calculations of the similar 
SUB(2,2) scheme not only yield similar results at J = 1/2 but the whole spectrum 
as a function of J [22].) These results of the triplet spectrum from such a low-order 
approximation seem to agree with a more substantial calculation at J = 1/2 by 
Shastry and Sutherland [23] . They obtained the spectrum of a soliton-like excitation 
with a minimum gap of 0.25 at g = and tt and a maximum gap of about 1 at 
q — 7r/2. Tonegawa and Harada's numerical calculations [20] confirmed the nonzero 
gap at J = 1/2 and in the nearby region. They predicted that the gap collapses at 
Jo ~ 0.3, while Haldane [24], who used a fermion representation, predicted this value 
to be about 1/6. Recently, Jq has been estimated by the conformal field theory to be 
about 0.2411 [25]. In any case, this gapless point may correspond to a phase transition 
from the dimerized phase to a critical phase similar to the Heisenberg model at J = 0. 
A more intriguing situation occurs for J > 1/2, where the triplet spectrum of 
both spin- wave theory and the SUB (2,1) scheme has a minimum at g = 7r/2. This 
certainly reminds us of the magneto-roton excitations in the fractional quantum Hall 
effects [26]. As J increases, the minimum (spin-roton gap) decreases and finally at 

(2) 

J = Jc , it collapses at q = 7r/2. Whether or not this suggests a phase change in the 
spatial periodicity of the system from double to four-fold, for example, is still unclear. 
The numerical calculations of the spin-spin correlation function [20] certainly show a 
more complicated feature for J > 1/2. In particular, the short-range four- fold Neel 
order (ttiittii ' ' is observed for J > 1/2, contrast to the case of J < 1/2 where 
the ordinary short-range two-fold Neel order (titi • • •) is observed [27]. As we know, 
at J = oo, the model Hamiltonian of Eq. (3.1) becomes two uncoupled Heisenberg 
chain with a double lattice spacing showing a four-fold spatial periodicity. Clearly, 
higher-order calculations are needed to obtain a clearer picture. 

To conclude this section, it should be pointed out that the dimerization order 
parameter, defined hy D = (si_i • Si) — (si • s^+i), can be easily obtained within spin- 
wave theory and the SUB2 scheme of the CCM. In spin-wave theory, for example, D is 
found to be nonzero in the region of < J < 1 and gradually diminish when J moves 
toward the two terminating points; but at the two terminating points {J — 0, 1), D 
diverges to — oo, implying a breakdown of spin- wave theory. The SUB2 scheme of 
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Figure 3. Schematic plots of the triplet excitation spectrum in the SUB(2,1) 
scheme at several values of J. 



the CCM, however, yields converging results even at the two terminating points, as 
it is the case in our previous CCM analysis for the ID anisotropic-Heisenberg model 
[18]. Since it also involves the ground bra state which is not manifestly the Hermitian 
conjugate of the ground ket state in the CCM, our CCM analysis for the dimerization 
order parameter and correlation functions will appear elsewhere. 



4. OTHER SYSTEMS 

4.1 The Spin-1 Chains 

Recently, the ID Heisenberg-biquadratic spin-1 chain has attracted much at- 
tention because it provides very rich and interesting quantum phases. The model 
Hamiltonian is given by 

H = cosOy^Si ■ Si+i + sinOy^isi ■ Si^iY , s = l (4-1) 



where the coupling between spins is parametrized by 9. The FM phase is restricted 
to the region of 7r/2 < 9 < 5tt/A, and the rest is non-FM. 

There are a number of exact results available at several values of 9. In particular, 
a,t 6 = — 7r/2, the system is exactly known to be dimerized with a nonzero gap and 
the corresponding order parameter is about 42% of the perfect dimer state [11,12]. 
At tanS = 1/3, the ground state is given by a homogeneous VB configuration with 
a nonzero energy gap but with no lattice symmetry breaking [7]. At 6 = 7r/4, the 
model is again integrable, the ground state clearly shows a triple spatial periodicity 
and the excitation spectrum becomes gapless at the lattice momentum q ~ and 2tt/3 
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Figure 4. Expectation values of the ID spin-1 Hamiltonian as a function of 
9 with respect to the three simple VB states: homogeneous (dotted), dimer 
(short dashed), and trimcr (long dash). Also shown are the results from 
finite-size exact calculations (solid). 

[28]. Based on these exact results, one tends to conclude that the system may show 
different phases representing by the three VB states respectively, namely the dimer 
state as shown in Fig. 1 (b), the trimer state in Fig. 1 (c) and the homogeneous VB 
state discussed in Ref. [7] . The expectation values of the Hamiltonian with respect 
to these three trial wavefunctions can be straightforwardly obtained as 

— = <( -cos6'+ fsine, dimer; (4.2) 

These values are shown in Fig. 4 as a function of 0, together with the numerical 
results from finite-size exact calculations [18]. One sees that the dimer state has 
lower energy than that of the homogeneous VB state for 9 < tan^^(— 1/2) k, —26.6°; 
for larger 9, however, the homogeneous VB state has lower energy. In particular, 
the homogeneous VB state is the exact ground state at tan 61 = 1/3 [7]. At even 
larger 9, it is interesting to see that the trimer state has the lowest energy. This 
occurs when 9 > tan~^(3/4) ?« 36.9°. The lower envelope of the three curves is in 
general quite close to the 'exact' results over the entire non-FM region. This crude 
approximation certainly seems to give a clear picture for the three-phase diagram of 
the spin-1 system, so far as the ground-state energy is concerned. Of course, the 
precise locations of the boundaries between these phases given here are not to be 
trusted because of the gross simplification. 
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From the above analysis, it is clear that one can extend our previous calculations 
to study the dimerization of the chain around the region oi 9 — — 7r/2 and to study 
possible trimerization about > 7r/4. Chubukov [29] applied a dinierized spin-wave 
theory using the Holstein-Primakoff bosonization to the Hamiltonian of Eq. (4.1) and 
indeed he found that over an extended region, the dinierized spin-wave excitations 
are stable. One certainly desires to obtain also other physical quantities, such as the 
ground-state energy, dimerization order parameter and the corresponding four-spin 
correlation functions, etc. The COM analysis described in Sec. 3 for the spin-i model 
can certainly provide a systematic means to obtain these physical quantities. 

The possible trimerization of the spin-1 chain was discussed in Ref. [13] where 
the equations of motion were derived for the creation and destruction operators AJ^q 
and Aq^ (n ^ 1,2,..., 26) with r denoting each of the trimcrs in Fig. 1 (c). After a 
truncation in the Hilbert space, namely restricting to n = 1, 2, . . . , 9, the trimerized 
spin-wave spectra were obtained. The lowest mode shows a nonzero gap associated 
with the trimerization, and this gap collapses at precisely 9 — -k /A and — 7r/2. In 
particular, a,t ~ 7r/4, the spectrum becomes gapless at lattice momentum q = Q and 
27r/3 with a spin-wave velocity of 3/\/5 w 1.342. This spectrum compares well with 
the exact result of Sutherland [28], which has a spin-wave velocity of \/2 7r/3 « 1.481. 
Kt 9 — 7r/2, where the system is known to become FM, a constant zero spectrum was 
obtained. Again, the CCM analysis using AJ^q and A^n {n — 1,2,..., 26) with the 
trimer model state should provide more systematic and reliable results. 

4.2 The Spin-i Frustrated Model on the Square Lattice 

The 2D spin-i frustrated Heisenberg model on the square lattice is described by 
the following Hamiltonian 

H ^-'^s,-s,+p + -J^Si-s,+p,, s = - (4.3) 

where i runs over all lattice sites, and p and p' over all nearest-neighbour and next- 
nearest-neighbour (diagonal) sites respectively. Because its possible relevance to the 
high-temperature superconductors, a variety of techniques has been applied to this 
model [14]. One now generally believes that the system shows the classical Neel- 
like order with the ordering wavevector Q = (tt, tt) for small J and the collinear 
magnetic order [Q = (0,7r)] for J « 0.65 or larger. Between these two phases (i.e., 
0.35 < J < 0.6), however, no magnetic order is observed. Although there is no clear 
consensus on the zero-temperature structure for this nonmagnetic region, the column 
dimerized phase shown in Fig. 5 has been proposed [14] . In particular, Chubukov again 
applied his dimerized spin- wave theory to the Hamiltonian of Eq. (4.3). His results 
seem to agree with the numerical calculations which suggest that the column dimer 
VB state may be stable around J = 1/2. But it is fair to say that a more systematic 
approach is needed before one can reach a definite conclusion on the dimerization of 
the 2D square lattice. It is quite straightforward to extend our CCM analysis for the 
ID spin-i frustrated chain described in Sec. 3 to the present 2D case. We will report 
these results soon. 
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Fig. 5. A column dimcr VB state for the 2D frustrated Heisenberg model 
on the square lattice. 



4.3 The Spin-1 Models on the Kagome Lattice 

Spin models on the kagome lattice are another group of frustrated systems be- 
cause the ground state of the classical Ising model on the kagome lattice has infinite 
degeneracy. In addition to their intrinsic theoretical interest, some spin models on 
the kagome lattice may have been realized in experiments. For example, in a layered 
compound Sr-Cr-Ga-0, the s = 3/2 Cr"^"*" ions form a stack of dense kagome lattices 
separated by more dilute triangular lattices [30]. And the spin-i Heisenberg model on 
the kagom,e lattice has been proposed to explain the interesting phenomena observed 
in the experiments with ^He atoms deposited on the graphite substrate [15]. 

Here the case of the spin-1 models on the kagome lattice is considered. It is useful 
to study the following trimerized Hamiltonian, 



H{J)= |E+^EH^-«^' ^ = 1' (4.4) 

where, as shown in Fig. 6, {ij) denote the solid bonds and {ij) the dashed bond. The 
symmetric model is given by J = 1. At J = 0, the perfect trimer VB configuration 
(i.e., sohd bonds in Fig. 6) is the exact ground state. Therefore one expects that 
the trimerized spin- wave theory discussed in Sec. 4.1 for the spin-1 chain should be a 
good approximation for at least small values of J. The method of equation-of-motion 
has been applied [31] to the Hamiltonian of Eq. (4.4) for the operators AJ^g and A^^ 
with n = 1,2,.. .,9, restricting the Hilbert space of each trimer (denoted by the 
new lattice vector r) to the first ten states. The trimerized spin-wave spectra have 
been obtained as functions of J. Unfortunately, the spectra are found to be stable 
only when J < 1/2, and the symmetric point J = 1 seems to be beyond this simple 
spin-wave approximation. 

However, similar to the ID spin-1 case discussed earlier, one can in general con- 
sider the spin-1 Heisenberg-biquadratic model on the kagome lattice. This model is 
given by adding a quadratic term to Eq. (4.4), 

iJ' = J'^(s,.s,+p)2, (4.5) 

where, as before, p denotes nearest-neighbour sites on the kagome lattice. From the 
experience of the ID spin-1 chain, one expects that the quadratic term of of Eq. (4.5) 
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Fig. 6. The spin-1 kagome lattice. 



may stabilize the trimer VB state over an extended region of J' and even at the 
symmetric point of J = 1. This work is in progress. 



5. CONCLUSION 

In this article, a microscopic approach to the quantum spin systems with an 
anticipated VBL long-range order is described. The perfect VBL consisting of in- 
dependent simple VBs is taken as the model state and the corresponding composite 
operators first developed by Parkinson are employed. Two approximation schemes 
are developed, firstly a spin-wave theory via bosonization transformation and sec- 
ondly the more systematic analysis within the framework of the CCM. The general 
formalism for the quantum correlations in the ground and excited states of the VBL 
systems are given. In particular, the simple ID spin-i frustrated model have been 
investigated in detail as a demonstration. The extensions of our approach to the 
spin-1 Heisenberg-biquadratic chain and to the 2D frustrated models on the square 
lattice and kagome lattice are also discussed. The preliminary results presented in 
this article are quite promising indeed. There is much more work to do. we wish to 
report our new results in the near future. 
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APPENDIX 
SYMMETRY THEOREM OF THE COUPLED-CLUSTER METHOD 

In Sec. 3.2, the SUB2 scheme of the CCM provides a solution which preserves 
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the rotational symmetry (i.e., Stotai = 0) of the model state \D), although one started 
with the operators in the sector of s^^^^^ — 0. We notice that this is true in a 
general application of the CCM. A brief proof of the following theorem is given in 
this Appendix: the CCM equations at any level of approximations always provide at 
least a solution which guarantees the symmetry of the model state if this symmetry 
is one of those belonging to the Hamiltonian. 

Let A be a symmetry operator, associated with which the model state |<J>) has 
an eigenvalue Aq, namely A|$) = Ao|$). Let Hamiltonian H commutes with A, 
[A, iif] = 0. Therefore, one has 

Aff|$) = Aoif|$>. (A.l) 

(I). Eigen operator representation. Let C^ and Cj" be the multi-configurational 
creation operators with respect to |$), with the set- indices / and J respectively label- 
ing the general multi-particle cluster configurations. The corresponding destruction 
operators are denoted as C/ and Cj respectively. Assuming that C^ commutes with 

A, 

[A, C+] = 0, V/, (A.2) 

but Cf does not. Instead, C^ has the following commutations, 

[A,C+]-A'i(J)C+, A'i(J)^0, VJ. (A3) 

Therefore one has for any positive integer n 

Aic+rm^xoic+rm, v/, (aa) 

and 

AC+|$)=Ai(J)C+|$), Ai(J)eeAo + A'i(J). (A5) 

The CCM correlation operator S is defined as 

5^^'5zC+ + ^5jC+, (A6) 

/ J 

where the prime implies that the identity term, Cq = 1, is excluded. The correlation 
coefficients {Sj, Sj} of Eq. (A. 6) are determined by the following sets of the coupled 
equations, 

($|C7e-^i?e^|$) = 0, V/ {^ 0), (A7a) 

($|Cje-^i?e^|$) = 0, VJ, {A.7b) 

where the expansion of the similarity-transformed Hamiltonian 

e-'^He^ = H+[H,S] + ^ [[H, 5] , 5] + • • • , (A.S) 

will terminate at a finite order in S* if iJ contains a finite number of destruction 
operators. 
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We notice that the states iJ|*), [ff, C/]|$), [[iJ, C/],C]^]|$), . . ., all have the 
eigen value Ao for A by Eqs. (A.l) and (A. 4). However, from Eq. (A. 5), the state 
Cj 1$) has a different eigenvalue, Ai(J) (^ Aq). Therefore they must be orthogonal 
to one another, namely 

($|Cji7|$) = 0, VJ (A9a) 

($|C.7[iJ,C/]|$) ==0, VJ,/ (AM) 

{^\Cj[[H,C+],C+m=0, yj,I,I' (A9c) 

etc. Using Eqs. (A. 8) and (A. 9), one immediately concludes that Eqs. (A. 7) have at 
least a solution given by 5j = for all J. If this solution is taken, the correlation 
operator S of Eq. (A. 6) preserves the symmetry of 1$). 

(II). Non- eigen operator representation, let Cy = (C/(0), C/(l), ..., Cy (n/ — 1)) 
be the multi-configurational creation operators with respect to 1$), with nj the di- 
mension of the symmetry A within a set /. Assuming Cj{n) do not have the com- 
mutations as in Eq. (A. 2) or Eq. (A. 3). Let the corresponding correlation coefficients 
be Si = (5/(0), 5/(1), ...,Si{nj — 1)). The CCM correlation operator is then written 
as 

S = Y! SrC+ = Y![Si{0)C+{0)+Si{1)C+{1) + - ■ ■+Si{ni-l)C+{ni-l)], (AlO) 

where the primes again imply exclusion of the identity term, / = 0. 

As before, the correlation coefficients {Si{i),i = 0, 1, ...,n/ — 1} are determined 
by the following set of the coupled equations 

($|C/(i)e^^iJe^|$) =0, i = 0, l,...,n/- 1 andV/ (^0). (All) 

Let B+ [= (B+(0),B/(l),...,B/(n/ - 1))] = T/ • C/, where T/ is a c-number 
(n X n) matrix and is chosen such that Bj{0) commutes with the symmetry operator 
A, namely 

[A,i?+(0)]=0, V/, (A12) 

but 

[A, B+{j)] = \]{I)B+{]), A;(/) ^0, J = 1, ..., n/ - 1 and V/. (A13) 

Multiplying Eq. (A. 11) by the corresponding elements of the c-number Hermitian 
matrix of T/, and after a simple summation, one derives the following equivalent 
equations 

($|B/(j)c-^i/e^|$) = 0, j==0, l,...,n/-l andV/ (t^O), (A14) 

where S can be equivalently written as 

S^S'i- B+, S'j = Si- Tj-\ (A.15) 

According to (I), there is at least a solution to Eq. (A. 14) hence to Eq. (A. 11) in 
which Sj{j) = for all j ^ 0. If this solution is chosen, the symmetry of the model 
state is preserved. Q.E.D. 
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